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Abstract 

, We investigate an electron in the plane interacting with the magnetic field due 

\Q ' to an electric current forming a localized rotationally symmetric vortex. We 

show that independently of the vortex profile an electron with spin antiparallel 
to the magnetic field can be trapped if the vortex current is strong enough. 
In addition, the electron scattering on the vortex exhibits resonances for any 
spin orientation. On the other hand, in distinction to models with a localized 
flux tube the present situation exhibits no bound states for weak vortices. 

P. 

Interaction of charged particles with a localized magnetic field has been a subject 
of interest for a long time, both from the theoretical and experimental point of 
view — see, e.g., ||GBG| , ptf and references therein. Such a field can have different 
sources, for instance, it may be induced by an electric current having one or more 
vortices. A lot of attention was payed in the last decade to vortex bound states in 
^ superconductors whose dynamics is governed by the Bogoliubov-de Gennes equation 

- cf. [|HKD| , gHJJSj, |G5| and a bibliography given in PTM[ . 



Another, much simpler example involves a Pauli electron interacting with a flux 
tube modelling a vortex magnetic field — it is is appealing, in particular, since it has 
been observed that vortices appear often in the probability current associated with 



mesoscopic transport — see [[ESSF|1 and the literature there. The system of a tube 



and an electron has been investigated in a fresh paper by Cavalcanti et al. [pFC 



who, however, seem to be unaware of another recent studies of the problem — cf. 
Mo|| and references therein. In these papers, the field is assumed to be constant 



within a circle and zero otherwise — the conclusion is then that in one spin state 
the electron can be always trapped by the vortex, independently of the magnetic 
flux value, as long as the effective gyromagnetic factor g* > 2 . This covers the 
physically important case of a free electron with g* = 2.0023. 

However, this claim depends substantially on the used magnetic field Ansatz. To 
illustrate this point we analyze in this letter the situation where the vortex current 
distribution represents the input. On one hand we are able to generalize the result 
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of [pVloj , |CFQ| by showing that vortices can trap electrons independently of their 



profile, and what is physically equally important, that they cause resonances in the 
electron-vortex scattering which get sharper as the vortex strength grows. On the 
other hand, the trapping needs in this more realistic setting certain minimum vortex 
strength (measured by the total circulating current or, say, the dipole moment of 
the field) to occur. To make things simpler we suppose that the vortex is centrally 
symmetric. We are convinced that the symmetry is not vital for the conclusions, 
but its absence makes the treatment technically much more complicated and we 
postpone the discussion of the general case to another paper. 

The dynamics of a nonrelativistic electron in the plane exposed to a perpendic- 
ular but nonhomogeneous magnetic field is given by the respective Pauli equation, 
i.e., by the Hamiltonian 

H = -L (-ihV + -A{x)) 2 + l 9f i B a 3 B{x) , (1) 



2m* V c / 2 

where m* is the effective electron mass, g the gyromagnetic factor, [ib = eh/2m e c 
the Bohr magneton, and the sign choice corresponds to the negative charge — e of 
the electron. In the rational units, 2m* — h — c — 1 this becomes 

H = (-zV + eA(x)) 2 + ~g*ea 3 B(x) (2) 

with the effective gyromagnetic factor g* = g — . Since H is matrix-diagonal, 
each spin state can be treated separately. Moreover, we shall consider the situation 
when the magnetic field is rotationally symmetric. This allows us to perform the 
partial-wave decomposition and to replace H by the family of operators 

*f } = - £i - H + > V t\r) ■= [eA{r) + ^ ± \g*eB{r) (3) 

on L 2 (M + , r dr) . To determine their spectral properties, we have to specify the 
involved functions, the angular component A(r) of the vector potential and the 
related magnetic field, B(r) = A'(r) + r~ 1 A(r) . 

We have said that they correspond a circulating electric current in the plane, 
which is supposed to be anticlockwise and having the angular component only, 
J(x) = \5(z) Jir)^ . Here r, ip are the polar coordinates in the plane; the total 
current is A J °° J(r) dr . The positive parameter A can be, of course, absorbed into 
J ; we introduce it as a tool to control the vortex "strength" . The current density 
is supposed to obey the following modest requirements: 

(i) J is C 2 smooth and non-negative, J(r) > 0, 

(ii) at the origin J(r) = ar 2 + 0(r 3 ) , 

(Hi) at large distances J(r) = 0(r~ 3 ~ e ) for some e > . 
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The vector potential in the plane z = is then also anticlockwise; its magnitude 



A{r) = A / dr'r'J{r') 



2tt 



cos ip' dip' 



(r 2 + r' 2 — 2rr' cos ip' 



,V2 



(4) 



is obtained by summing the contribution from all the circular current lines |Ja 

4r (2 - p 2 )K(p 2 



where p : 



A(r) = A / J if) 



4rr' 
(r+r 1 ) 2 



mp 2 



dr' 



(5) 



and K, E are the full elliptic integrals of the first and the second 



kind, respectively. Using ||AS| , 17.3.29,30] (pay attention to a misprint there), one 
can also cast (IST) into the form 



A(r) 



4A /„ J ^T < [ K w- E wr- 



(6) 



where we have used the usual shorthands, r 



minfr, r'l and r- 



max r, r 



the same result follows directly from (^) and |pR| , 3.674.3]. In view of (%) the integral 
is finite for every r , because E(() is regular at £ = 1 and -f^(C) has a logarithmic 
singularity there. 

Let us denote the Pauli Hamiltonian (H) with the vector potential (^) by if (A) ; 
the symbol H\ (A) will be used for its spin and orbital momentum components. 
Our main result is then the following: 

Theorem. Under the stated assumptions, a(H(X)) = [0, oo) for |A| small enough. 
On the other hand, Hq \\) has a negative eigenvalue for a sufficiently large A . 

Proof: The first claim has to be checked for any (X) . By (i), the potentials 
are C l smooth; they decay at infinity as 



K (±) M 



P 2 If ~x. n* 

1 + Xem^- + 0(r" 3 -) 

ny £ jiO 



(7) 



where m := tt / °° J(r') r' 2 dr' is the dipole moment of the current for A = 1 - - cf. 
fja|. Consequently, cr e < jS (if^ ± ' ) (A)) = [0, oo) following [[R5| , Sec. XIII. 4]. We rewrite 
the potentials into the form 



V^\r) = [XeA 1 (r) + 



(8) 



where the the indexed magnetic field refers to the value A = 1 . Since H[ ; (A) is 
just the s-wave part of the two-dimensional Schrodinger operator with the centrally 
symmetric potential (§), it is sufficient to check that the latter has no negative eigen- 
values. If C = the potential decay allows us to apply Thm. 3.4 of |S| by which 
a negative eigenvalue exists for small positive A if and only if J °° (r) r dr < . 
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However, the flux through the circle of radius r is 2tt f Q B(r')r'dr' = 2irrA(r) , so 
the second term in (||) does not contribute and the integral is determined by the first 
one which is positive for any nonzero A . If £ 7^ the decay is too slow, but this diffi- 
culty is easily overcome. We replace the first term, e.g., by (\eAi(r) + £\ 6(r — r) 
with a positive r and obtain the absence of a negative eigenvalue for a small A ; 
the same is by the minimax principle true for the original operator. 

For the existence claim the behaviour of the potential around the origin is vital. 
We shall write the vector potential in the form 

A{r) = A/ir + a (r) , n := / J(r') — . (9) 

Jo r 

The behaviour of «o follows from (Q) and the following estimate on the difference 
of the elliptic integrals, 

Y < K (0 - E (0 < \ (l - I ln(l - C) + b(A (10) 



for a sufficiently large b > , which is a straightforward consequence of |AS 
17.3.11,12] and [^, 8.123.2]. On the lower side we thus get 

7r f r r r dr 1 3 

otoir) > -it / J(r')r' 2 dr' — irr / J(r') — > nar 2 + 0(r 3 ) . 

r l Jo Jo r 4 

The upper bound is similar. In view of (i)-(iii), J(r) < dr for a suitable a ; in the 
logarithmic term we employ the Taylor expansion which gives 

37m r°° ~ r 2 i , _ 3™ 2 ^ 1 

Together we find ao(r) = 0(r 2 ) . This further implies 

B(r) = 2Xfi + Po(r) , (3 (r) := a' (r) + l -a Q {r) e 0(r) . (11) 

The operator Hq \\) can be therefore written as 

^ _) (A) = - ~ + (Ae^r + aoW)) 2 - \ Xeg*(2^ + (3 (r)) . (12) 

Using the rescaled variable u := ry/X we find it to be unitarily equivalent to the 
operator 

\A\ with A x = A + W\ (13) 

on L 2 (M + ,udu) , where 

Ao - ~ -4- ~ 9*e^ + e 2 n 2 u 2 (14) 
du A u du 



and 



W x (u) 



2%/Ar 2 



e /juao 




+ Ae 2 al 





(15) 



We have clearly a ess (A\) = <r ess (\A\) = [0, oo) for any A > , while A Q as the s- 
wave part of the two-dimensional harmonic oscillator has a purely discrete spectrum. 
Despite the fact that the perturbation is large (the maximum of W\ grows linearly 
with A ) , one may attempt to employ the asymptotic perturbation theory. 



Unfortunately, Thms. VIII. 3. 11, 13 of |KaH cannot be applied directly, because 
the family {W\} is not monotonous. Instead we use the fact that W\ tends to zero 
pointwise as A — > oo , since the above estimates yield 



\W x {u)\ < (2e 2 fxc a u 3 + ^g*ec p v\\- 112 + 



e 2 ciu A \- 1 



(16) 



for some positive c a , cp . The family {A\} can be estimated from below because 
the potentials have in view of flT5D a uniform lower bound. Hence one can choose 
£o < which belongs to p(A\) for all A , and the resolvents form a bounded family 
of positive operators, (A\ — £) _1 < (£o — _1 > f° r an Y £ < £o • Next we use a 
trick based on the resolvent identity: for a vector / = (Aq — with a fixed 
g E C^(R + ) we have 

\(A X - f - (Ao ~ f\\ = I^A-fl" 1 ^! < (to- 0' 1 \\W x g\\ - (17) 

as A — >• oo in view of ( pp and the compact support of g . However, the family of 
such / is dense in L 2 (M + ,udu) , so A\ — > A in the strong resolvent sense. 

This allows us to employ Thm. VIII. 1.14 of | |Ka| by which to any v n e a p (A Q ) 
there is a family of u n (X) G o-(A\) such that u n (\) — > v n as A — > oo . It is therefore 
decisive that the unperturbed eigenvalue is stable in the sense of [|Ka , which means 
negative in our case. The spectrum of Aq is given explicitly by 



e\i (An + 2 - 



0,1, 



(18) 



so the condition is satisfied for n = if g* > 2 (as in [|UF(J|| , the next eigenvalue 
comes to play for g* > 6 which is physically not appealing). □ 

We finish the letter by remarks on extensions of the result and related topics: 
(a) The difference in the weak-coupling behaviour comparing to |Mo, CFC 



is not 



surprising. In their case / °° Vq(t) rdr is dominated for small A > by the 
negative term due to the well. In reality, however, the magnetic field flux lines 
are closed in M 3 , so the well is compensated by a repulsive tail, small but 
extended, which prevents the trapping. 

(b) The asymptotic perturbation theory yields also the ground state behaviour as 
A — > oo . By [|Ka , Thm. VIII. 2. 6] the leading-order correction to Uq is 




a'^)\- l ' 2 + 0{\- 
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where ipo is the ground-state eigenfuction of the two-dimensional harmonic 
oscillator. However, ad'(O) = , so the ground state of the original operator 
H^~\\) behaves as -Ae/i(#*- 2) + 0(1) . 

(c) Large A give rise to an orbital series of bound states: the above also argument 

works for B.\ (A) with i = —1,-2,.... The potential in ( |I4D is replaced 
at that by e 2 fi 2 u 2 + £ 2 r~ 2 + e\i{2£ — g*) , and one looks for negative eigenval- 
ues among u Hj i = e/i (An + 2(\£\ +£) + 2 — g*) . The critical A at which the 
eigenvalue emerges from the continuum is naturally ^-dependent. 

(d) Positive eigenvalues of A are unstable in the sense that they disappear in the 

continuum once the perturbation ( ]T5|) is turned on. Following |[Ka| , Sec. VIII. 5], 
however, they give rise to spectral concentration as A — > oo which is mani- 
fested by resonances in electron scattering on the vortex. Knowing the shape of 
the potential barrier, one can compute their widths which vanish exponentially 
fast with A . For a fixed I , the number of resonances grows asymptotically lin- 
early with A , because the eigenvalues of the operator (|T4]) are equally spaced 
and the top of the potential barrier in A\ is asymptotically linear in A . Due 
to the same reason, resonances exist at large A for both spin signs. Notice also 
that the existence of resonances is not restricted by the value of the gyromag- 
netic factor, and therefore they may be observed in semiconductor systems 
where \g*\ is typically less than one, in some cases even of order of 10 -2 . 

(e) The above mentioned resonance scattering is a purely quantum effect. A clas- 

sical electron can be, of course, trapped in the current vortex if it is placed 
inside the potential barrier (|3|) with the energy less than its top (since there is 
no spin in this case, g* = , the well bottom is at zero). On the other hand, 
an electron of the same energy outside the barrier gets scattered without the 
possibility of entering temporarily the interior of the vortex. 

(f) In addition to smooth current distributions discussed above, the case J(r) = 

5(r — R) is of a practical interest. As an illustration, imagine two adjacent 
thin films, one supporting a free electron gas while the other is equipped 
with a mesoscopic ring in which a persistent current circulates. The effective 
potential is now given explicitly by (Q) and @; in distinction to the smooth 
case it has a singularity of the type (r— i?) -1 and the operators if^(A) need 
a regularization. Since the 5 function above is an idealization of a sharply 
localized distribution, it is reasonable to choose for this purpose from the 
family of all admissible procedures a scheme based on the principal value of 
the singular potential ||N Z] , |Ku|| . 

(g) The reader may wonder whether the preferred spin orientation of the ground 
state does not conflict with the second theorem of Aharonov and Casher ||AC|| . 
The present model offers a good illustration that this claim should be taken 
cum grano salis indeed (in distinction to the first one which is sound — cf. 
||Th|| ). Recall that the Pauli Hamiltonian (0) can be factorized into the product 
(III- io- 3 U 2 )(Jl 1 + K7 3 n 2 ) , where II := -iV +eA; if * solves if* = Ety 
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then \1/ := (IIi+ ias^)^ solves the equation with spin term sign switched. 
The nonzero component of the rotationally symmetric spinor \l/ turns at that 
into — ie l<p (ijj' (r) — eA(r)ip(r)) . For the ground-state eigenfunction we have 
ip{0) 7^ while i[>'(0) = due to the smoothess of \& . However, the vector 
potential satisfies (0) with a positive coefficient. To belong to the domain of 
H , the nonzero component of \1/ must be smooth at the origin. This requires 
^"((^(O) -1 = Ae/i but the Ihs is negative; recall that ip approaches the 
harmonic-oscillator ground state for large A . 

(h) Notice finally another difference. The model with a flux tube can have a natural 
"squeezing limit" in terms of Aharonov-Bohm Hamiltonians with a pointlike 
magnetic flux ]Mc| |DS| because the vector potentials coincide outside the 



flux tube and the attractive part tends to a 5-well (with the exception of a 
single one, the resonances mentioned in (d) are lost at that as it is usual in 
such situations — cf. [ |AGHH| |). On the other hand, the present case is more 
complicated being essentially three-dimensional as far as the magnetic field is 
concerned. 

Acknowledgment 

P.E. thanks for the hospitality extended to him at C.P.T. where this work was done. 
The research has been partially supported by GAAS under the contract A1048801. 



References 

[AS] M.S. Abramowitz, I. A. Stegun, eds.: Handbook of Mathematical Functions, 

Dover, New York 1965. 
[AT] R. Adami, A. Teta: On the Aharonov-Bohm Hamiltonian, mp_arc 97-37. 
[AC] Y. Aharonov, A. Casher: Ground state of a spin-1/2 charged particle in a 

two-dimensional magnetic field, Phys. Rev. A19 (1979), 2641-2642. 
[AGHH] S. Albeverio, F. Gesztesy, R. H0egh-Krohn, H. Holden: Solvable models in 

quantum mechanics, Springer, Berlin 1988. 
[CFC] R.M. Cavalcanti, E.S. Fraga, C.A.A. de Carvalho: Electron localization by 

a magnetic vortex, Phys. Rev. B56 (1997), 9243-9246. 
[DS] L. Dabrowski, P. Stovicek: Aharonov-Bohm effect with the 5 type inter- 
action, J. Math. Phys. 39 (1998), 47-62. 
[ESSF] P. Exner, P. Seba, A.F. Sadreev, P. Stfeda, P. Feher: Strength of topologi- 

cally induced magnetic moments in a quantum device, Phys. Rev. Lett. 80 

(1998), 1710-1713. 

[GBG] A. Geim, S.J. Bending, I.V. Grigorieva: Asymmetric scattering and diffrac- 
tion of two-dimensional electron at quantized tubes of magnetic flux, Phys. 
Rev. Lett. 69 (1992), 2252-2255. 

[GR] I.S. Gradshtein, I.M. Ryzhik: Tables of Integral, Series, and Products, Aca- 
demic Press, New York 1980. 



7 



[GS] F. Gygi, M. Schluter: Angular band structure of a vortex line in a type-II 
superconductor, Phys. Rev. Lett. 65 (1990), 1820-1823. 

[HIM] N. Hayashi, M. Ichioka, K. Machida: Effects of gas anisotropy upon the elec- 
tronic structure around a superconducting vortex, Phys. Rev. B56 (1997), 
9052-9063. 

[HRD] H.F. Hess et al.: Scanning-tunneling-microscope observation of the Abri- 
kosov flux lattice and the density of states near and inside a fluxoid, Phys. 
Rev. Lett. 62 (1989), 214-216. 

[Ja] J.D. Jackson: Classical Electrodynamics, John Wiley, New York 1962. 

[Ka] T. Kato: Perturbation Theory for Linear Operators, Springer, Heidelberg 
1966. 

[Ku] P. Kurasov: Response to "Comment ..." by Fischer, Leschke and Miiller, 

J. Phys. A30 (1997), 5583-5589. 
[Mo] A. Moroz: Single-particle density of states, bound states, phase-shift flip, 

and a resonance in the presence of an Aharonov-Bohm potential, Phys. 

Rev. A53 (1996), 669-694. 
[NZ] H. Neidhardt, V.A. Zagrebnov: Towards the right Hamiltonian for singular 

perturbations via regularization and extension theory, Rev. Math. Phys. 8 

(1996), 715-740. 

[RS] M. Reed, B. Simon: Methods of Modern Mathematical Physics, IV. Analysis 

of Operators, Academic Press, New York 1978. 
[SHDS] J.D. Shore, Ming Huang, A.T. Dorsey, J. P. Sethna: Density of states in 

a vortex core and zero-bias tunneling peak, Phys. Rev. Lett. 62 (1989), 

3089-3092. 

[Si] B. Simon: The bound state of weakly coupled Schrodinger operators in one 

and two dimensions, Ann.Phys. 97 (1976), 279-288. 
[St] P. Stovicek: Scattering on a finite chain of vertices, Duke Math. J. 76 

(1994), 303-332. 
[Th] B. Thaller: The Dirac equation, Springer, Berlin 1992. 



8 



